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Abstract 


A  recently  proposed  equation  of  state  is  further  examined  for  compatibility  with  prior  theory. 
It  is  found,  by  adopting  a  more  realistic  functional  relationsh^  for  one  of  the  thermodynamic  state 
variables  and  eliminating  higher-order  terms  from  the  model,  that  three  historical  theories  for 
determination  of  the  Gruneisen  function  fall  directly  out  of  the  recently  proposed  model.  One  is 
led  to  conclude  that  the  newer  model  captures  the  essence  of  aU  of  the  older  theories,  each  of 
which  being  applicable  under  an  appropriate  circumstance. 
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1.  BACKGROUND 


The  estimation  and  measurement  of  the  Grtineisen  function  for  metals  is  an  important  but 
difficult  task.  Simple  formulas  exist  for  its  estimation,  in  terms  of  specific  heat  and 
compressibility,  or  alternately  in  terms  of  derivatives  of  the  cold  {i.e.,  zero  degree)  compression 
curve,  or  still  again  in  terms  of  derivatives  of  the  lattice  frequency.  Depending  on  the  validity 
of  various  assumptions  about  the  material,  some  of  these  definitions  may  be  redundant. 
However,  when  applied  to  actual  substances,  these  methods  can  produce  estimates  that  vary  by 
significant  amounts  (100%,  in  some  cases  [1]),  even  at  ambient  conditions. 

Some  variation  is  due  to  the  fact  that  the  Gruneisen  assumption,  which  has  the  Griineisen 
function  independent  of  temperature,  is  only  an  approximation.  It  is  well  known  [2,3],  for 
example,  that  at  low  temperatures  (at  fractions  of  the  Debye  temperature)  the  Gruneisen 
parameter  is  indeed  temperature  dependent.  Grodzka  [1]  presents  data  showing  a  significant 
thermal  variation  in  the  Gruneisen  parameter  at  higher  temperatures  as  well,  for  some  materials. 
Other  variations  in  the  estimate  arise  out  of  the  differences  in  modeling  of  the  propagation 
velocities  associated  with  lattice  vibration.  Three  widely  referenced  models  of  this  variety,  which 
provide  not  only  an  ambient  value  but  also  a  volume-dependent  estimate  of  the  Gruneisen 
function,  are  those  of  Slater  [4],  Dugdale  and  MacDonald  [5],  and  the  free- volume  theory  [6]. 
Pastine  [7]  and  Vashchenko  and  Zubarev  [6]  have  attempted  to  generalize  results  by  showing  that 
the  variation  of  Poisson’s  ratio  with  volume  can  account  for  some  of  the  differences  in  the 
various  models  of  the  Gruneisen  parameter.  Unfortunately,  experimental  knowledge  of  how 
Poisson’s  ratio  varies  with  compression  and/or  temperature  is  as  elusive  as  knowledge  of  the 
Gruneisen  function  itself. 

Historically,  the  estimation  of  the  Gruneisen  function  was  further  complicated  by  the 
paucity  of  cold-compression  data.  Available  high-pressure  data  were  almost  exclusively  that  of 
shock  transition,  and  so  the  cold-compression  curves,  needed  by  the  various  models  in  order  to 
estimate  the  Griineisen  function,  were  necessarily  back-extrapolated  from  shock  data  using  a 
Gruneisen  function  that  was  a  priori  unknown  {e.g.,  see  Walsh  et  al.  [8]). 
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To  some  extent,  this  uncertainty  in  the  functional  behavior  of  the  Gruneisen  function  has 
persisted  to  the  present.  For  example,  Segletes  [9]  examined  the  Gruneisen  formulations  present 
in  five  of  the  widely  used  hydrocodes,  which  were  developed  by  the  U.S.  Department  of 
Energy  (DOE)  national  labs  and  the  U.S.  Department  of  Defense  (DOD),  as  well  as  industry. 
In  examining  these  codes,  which  are  utilized  to  model  the  high-compression  deformations 
associated  with  high-velocity  impact  and/or  explosive  loading,  it  was  noted  that  a  common 
Gruneisen  formulation  was  not  employed.  Not  only  were  the  Gruneisen  formulations  different 
from  code  to  code,  but  Segletes  was  able  to  show,  in  some  cases,  that  thermodynamic  instabilities 
would  result  because  of  an  incompatible  selection  of  Gruneisen  and  Hugoniot  functions. 

In  a  significant  work,  which,  on  the  surface,  might  not  appear  to  shed  light  on  the 
problem  of  estimating  the  Gruneisen  function.  Rose  et  al.  [10]  proposed  (in  1984)  a  universal 
cold  curve,  which  was  able  to  provide  a  generic  function  for  the  lattice  binding  energy,  which 
could  be  appropriately  scaled  to  match  theoretical  and  experimental  data  for  many  metals.  Since 
the  work  of  Rose  and  his  colleagues,  others  have  continued  in  the  area  as  well,  notably  the  same 
team  led  by  Vinet  [11,12]  and,  more  recently,  a  group  led  by  Baonza  [13,14].  These  later  efforts 
succeeded  in  introducing  temperature  effects  and  thus  produced  a  Griineisen-style  equation  of 
state.  However,  in  all  of  these  more  recent  efforts  cited,  the  function  governing  the  reference- 
compression  curve  was  mutually  exclusive  of  the  function  used  to  model  the  thermal  effects.  In 
the  works  of  Vinet  et  al,  the  variation  of  the  bulk  modulus  was  fitted  for  a  number  of  materials, 
which  effectively  parlayed  into  a  Gruneisen  function.  In  the  case  of  Baonza  et  al,  a  particular, 
commonly  employed,  functional  form  of  the  Gruneisen  function  was  used — ^namely,  TIV  equals 
a  fitted  constant.  While  all  of  these  works  are  worthy  efforts  in  their  own  right,  their  primary 
contribution  lies  in  an  area  other  than  the  theory  of  the  Gruneisen  function. 

Segletes  [15,16]  recently  proposed  an  equation  of  state  for  metals,  which  captures,  in  the 
absence  of  phase  change,  the  behavior  of  metals  into  the  megabar  range.  The  model  was 
likewise  inspired  by  the  universal  cold-curve  work  of  Rose  et  al  [10].  However,  unlike  the 
efforts  of  Vinet  et  al  [1 1,12]  and  Baonza  et  al  [13,14],  in  which  the  function  governing  the  cold 
compression  was  mutually  exclusive  of  the  function  used  to  model  the  thermal  effects,  the  recent 
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work  of  Segletes,  by  expressing  the  binding-energy  function  in  terms  of  lattice  frequency,  as 
opposed  to  lattice  compression,  unified  the  functions  governing  cold-compression  and  thermal 
effects.  In  essence,  a  specification  of  the  cold-compression  behavior  in  Segletes’  model 
automatically  yields  the  Griineisen  function,  and  conversely,  a  specification  of  the  thermal  effects, 
by  way  of  the  Griineisen  function,  automatically  yields  a  cold-compression  curve. 


Segletes’  equation  is  of  the  Griineisen  variety  and  is  given  by 


pyif-E 


+  /«:(«■-  l)(0/©o)*^  ln(©/0o)} 


(1) 


where  p  and  E  are  the  pressure  and  specific  internal  energy,  respectively,  and  the  constant  E^  is 
the  specific  binding  energy  of  the  lattice.  The  terms  ©o  and  ©  are  the  ambient  and  current  values 
of  the  characteristic  temperature  of  the  lattice,  which,  according  to  the  tenets  of  Griineisen  theory, 
is  assumed  to  be  a  function  of  volume  only  and  is  proportional  to  the  maximum  lattice 
vibrational  frequency,  v.  The  concepts  of  characteristic  temperature  and  lattice  frequency  derive 
from  the  early  work  on  specific  heats  at  low  temperatures,  by  Einstein,  Debye,  Bom  and  von 
Karman,  and  others.  In  the  current  context,  the  characteristic  temperature  serves  as  a 
macroscopic  parameter  that  permits  the  vibrational  frequency  of  the  lattice  (a  microscopic 
quantity)  to  be  inferred — ^namely,  the  theory  of  specific  heats  has  these  two  parameters  directly 
proportional.  Returning  to  eqn.  (1),  the  parameter  \|r  is  the  thermodynamic  state  variable  equal 
to  the  ratio  of  specific  volume  to  the  Griineisen  parameter,  V/F,  introduced  by  Segletes  [17,18] 
for  convenience  in  manipulating  the  thermodynamic  equations.  The  constant  K  is  given  by 

K  =  _S_  ,  (2) 

To^/ir 


which  incidentally  is  l/CSFo)  times  the  anharmonicity  of  the  lattice.  Values  for  K  have  been 
observed  to  cluster  in  the  2/3  to  4/3  range  for  metals,  which  is  compatible  with  the  model’s 
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behavior  becoming  greatly  idealized  (i.e.,  linear  in  0/0o)  when  the  constant  K  takes  on  a  value 
of  unity.  Q  is  the  bulk  sound  speed  at  absolute  zero  temperature  and  pressure,  and  Fq  is  the 
ambient  value  of  the  Griineisen  parameter. 

In  this  model,  the  cold  curve  (i.e.,  zero-degree  isotherm)  is  given  by 


=  E,[l  -[1  -.S:in(0/0o)](0/0o)^} 


(3) 


and 


-  (0/0o)^  ln(0/0o)  ,  (4) 

'll/ 

where  and  E^  are  the  cold  pressure  and  specific  cold  energy,  respectively. 

In  macroscopic  terms,  the  Griineisen  function,  F,  is  defined  by  F  =  Wt|r  =  V(dpldE)y. 

From  microscopic  arguments,  however,  the  characteristic  temperature,  0  (alternately,  the 
maximum  lattice  vibrational  frequency,  v),  is  related  to  the  tir  variable  through  the  relationship 
that  defines  the  Griineisen  parameter: 


V  ^  0'  ^  -JL 

V  0  tir 


(5) 


where  primes  (here  and  throughout  the  paper)  denote  ordinary  differentiation  with  respect  to 
specific  volume.  In  the  general  case,  the  differentiation  in  eqn.  (5)  would  be  partial  at  constant 
temperature.  However,  for  Griineisen  materials,  in  which  F  =  F(V)  only,  the  derivative  in 
eqn.  (5)  becomes  ordinary. 
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It  is  interesting  to  note  the  distinction  between  Segletes’  equations  and  those  of  an  ideal, 
harmonic  oscillator,  which  is  often  used  to  approximate  a  crystal  lattice  for  purposes  of 
theoretical  derivation.  In  the  case  of  Planck’s  harmonic  oscillator  at  zero  temperature  [19], 
dEJdM  is  constant.  No  lattice  is  truly  harmonic,  of  course,  or  else  its  thermal  expansion 

coefficient  and  Griineisen  function  would  both  be  identically  zero.  From  eqns.  (3)  and  (5),  it  can 
be  seen,  for  the  idealized  case  where  K  is  unity,  that  the  current  model  would  have 
dEJdv  oc  InCv/Vg) .  A  theoretical  explanation  for  this  distinction  is  not  yet  postulated. 

Once  the  i);  function  is  defined,  integration  of  eqn.  (5)  gives  the  characteristic  temperature, 
which  is  the  primary  equation-of-state  variable  in  eqn.  (1).  In  the  original  papers, 
Segletes  [15,16]  incorporated  a  linear  relationship  for  \|;(V),  based  on  results  of  his  prior  work 
on  thermodynamic  stability  [9,17,18],  given  in  general  form  by 

=  i-(ry)(i-wy^  .  (6) 

The  ambient  values  of  specific  volume  and  Griineisen  parameter  determine  \|ro,  while  the 
constant  was  the  sole  material  parameter  that  was  actually  fitted  to  the  cold  curve  of 

each  metal.  The  fitted  values  of  (r(,\|/),  for  nine  metals,  ranged  from  0.43  to  0.88,  with  a  vast 

majority  clustered  in  the  vicinity  of  0.8.  In  terms  of  compression,  p,  equal  to  (Vo-V)IV,  the  linear 
translates  to  a  Griineisen  function  of  the  form  r=ro/(l+Pp),  with  P  as  the  fitting  parameter. 

This  assumed  linear  relationship  produced  an  excellent  match  to  empirical  data  and  first- 
principles  computations  over  a  very  wide  range  of  specific  volumes.  However,  for  extremely 
large  compressions,  this  linear  relationship  can  no  longer  be  expected  to  hold,  as  it  is  later  shown 
that  the  variable  must  vanish  in  the  high-compression  limit  (V  — >  0),  whereas  the  linear  form 

of  eqn.  (6)  will  not  do  so  for  values  of  (FoV)  other  than  unity. 
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In  this  report,  a  power  law  for  the  \|r(V)  parameter,  to  replace  the  linear  relationship 
originally  employed,  is  adopted.  Though  this  choice  of  function  still  suffers  deficiencies  at  very 
large  compressions,  a  number  of  positive  outcomes  result  in  the  context  of  Segletes’ 
model  [15,16]  with  this  change:  1)  the  \|f(V)  parameter  vanishes  with  diminishing  specific 
volume,  as  is  required;  2)  binding-energy  expansion  data  for  actual  metals  are  better  fit  with  the 
power  law,  as  compared  to  the  linear  assumption;  and  most  importantly,  3)  Segletes’  equation 
of  state  is  shown,  by  neglecting  higher-order  terms,  to  reduce  to  the  model  of  Slater  [4],  or 
Dugdale  and  MacDonald  [5],  or  to  the  ffee-volume  theory  [6],  depending  on  the  value  of  the 
exponent  used  in  the  power  law  for  \(r(V). 

2.  THE  HIGH-COMPRESSION  LIMIT  FOR  \|r(V) 

With  the  linear  relation  for  \|;(V),  given  by  eqn.  (6),  the  high-compression  limit  for  tir 
approaches  a  value  of  (1  -  FoV)?  which  would  typically  lie  in  the  0.2  •'V|rQ  range  for  metals. 

To  show  what  high-pressure  limit  is  appropriate  for  \|/‘,  in  the  context  of  the  new  model, 
theoretical  values  for  cold  pressure  and  energy  using  Thomas-Fermi-Dirac  (TFD)  theory  for 
several  metals  were  inserted  into  the  new  model  in  order  to  back  out  values  for  In  this 
manner,  the  high-pressure  limiting  functional  behavior  for  t|f,  which  would  be  required  to  match 
TFD  theory,  was  obtained.  Values  for  the  TFD  data  were  obtained  using  the  analytical  fitting 
equations  described  by  Kerley  [20]. 

It  is  not  asserted  here  that  the  new  equation  of  state  can  or  should  capture  the  high- 
compression  atomic  interactions  modeled  with  TFD  zero-temperature  theory.  Rather,  it  merely 
makes  sense,  firom  an  application  point  of  view,  that  the  current  model  should  blend  smoothly 
into  TFD  theory,  so  that  a  computational  implementation  can  smoothly  transition  between 
applicable  theories. 

Results  of  the  TFD  fits  for  \|;  are  shown  in  Figure  1,  for  several  metals.  Superimposed 
on  the  graph  are  lines  of  constant  which  characterized  the  original  linear  form  employed 
by  Segletes  [15,16].  A  finite  limiting  value  of  \|r  implies  that  F  approaches  zero  in  the  high- 
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Figure  1.  Limiting  11/“  values  for  several  metals  (solid  lines),  obtained  by  inserting  TFD  cold- 
curve  values  into  equation  of  state  of  Segletes.  Superimposed  (dashed  lines)  are  limiting  \|f 
values  associated  with  the  linear  \|r  assumption. 
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compression  limit,  which  is  known  not  to  be  the  case.  The  TFD  results  indicate  that  \|r  must  drop 
to  zero  with  decreasing  specific  volume,  which  indicates  that  the  limiting  value  of  T  is  nonzero, 
but  instead  related  to  the  limiting  slope.  Regardless  of  material,  the  TFD  results  are  roughly 
linear,  having  a  slope  of  approximately  2  in  the  very-low-volume  limit,  with  the  slope 
Himinighmg  slightly  with  increasing  volume.  A  finite  limiting  slope  of  value  m  in  Figure  1 
implies  a  limiting  value  of  the  Grfineisen  parameter  of  Fo/m.  The  fact  that  (in  the  case  of  m»2) 
ro/2  doesn’t  equal  1/2  or  2/3,  as  might  be  expected,  merely  indicates  that  the  current  model  fails 
to  capture  the  essence  of  TFD  theory,  which  is  no  great  surprise  since  the  model’s  equations  are 
not  those  of  TFD  theory.  In  any  event,  the  linear  treatment  of  \|r,  originally  employed,  does  not 
blend  smoothly  into  the  realm  where  TFD  theory  governs.  Another  functional  choice  for  \(r  might 
be  more  appropriate. 

3.  A  POWER  LAW  FOR  \ir(V) 

As  an  alternative  to  the  linear  relationship  for  \|r(V),  given  by  eqn.  (6),  a  power  law  of  the 

form 


=  (y/v,r  (7) 

is  investigated.  The  power-law  form  for  \|;(V)  retains  the  mathematical  simplicity  necessary  to 
manipulate  the  governing  equations  conveniently,  and  despite  a  few  deficiencies  that  will  be 
pointed  out,  it  seems  to  be  a  step  in  the  right  direction.  In  terms  of  compression,  p,  the  power- 
law  'V|/'  translates  to  a  Gruneisen  function  of  the  form  r=ro/(l+p)^ 

Figure  2  shows  the  same  TFD  results  as  presented  previously,  except  that  several  contours 
of  a  power-law  'Vl/‘(V)  are  superimposed  on  the  graph,  rather  than  the  linear  contours  of  Figure  1. 
At  first  glance,  the  power-law  \|/'(V)  might  seem  to  offer  only  minimal  benefit  over  the  linearized 
form.  Both  forms  cross  over,  rather  than  asymptote,  to  the  TFD  curves.  Also,  the  power-law 
curves,  though  reducing  to  a  value  of  zero,  as  is  required,  do  so  with  an  mfinite  limiting  slope, 
so  that  the  limiring  value  of  F  is  not  positive,  but  remains  zero.  On  the  positive  side,  however. 
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V/Vo 


Figure  2.  Limiting  \j/'  values  for  several  metals  (solid  lines),  obtained  by  inserting  TFD  cold- 
curve  values  into  equation  of  state  of  Segletes.  Superimposed  (dashed  lines)  are  limiting  \(/' 
values  associated  with  the  power-law  assumption. 
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the  curve  corresponding  to  an  exponent  of  0.8  does  blend  nicely  into  the  TFD  curves.  And  of 
the  nine  metals  examined  by  Segletes  [15,16],  seven  had  values  for  (corresponding  roughly 

to  exponent  x)  in  the  0.75  to  0.88  range,  where  the  power  law  fit  approximates  the  TFD 
asymptote  very  well. 

As  far  as  actual  cold-compression  data,  the  power  law  and  the  linear  assumption  for  \it(V) 
are  nearly  identical  for  values  of  0.8  over  relative  volumes  from  unity  down  to  0.5.  fri  this 

region,  the  power  law  provides  a  slightly  more  concave-upward  curvature  to  the  cold  curve, 
which  can  be  mitigated  by  choosing  a  value  of  exponent  x,  slightly  below  the  fitted  value  of 

Fotir'  (e.g.,  whereas  rQ\|r'=0.76  fits  stainless  steel,  a  value  of  x=0J2  fits  comparably  well).  Over 

the  range  of  available  cold-compression  data  [10,21,22],  the  power-law  and  linear  tirCV) 
formulations  can  be  made  virtually  indistinguishable  with  appropriate  selection  of  x,  both 
resulting  in  excellent  fits  to  the  data. 

On  the  expansion  end  of  the  curve,  the  power-law  \|r(V)  is  seen  to  provide  an  improved 
correlation  to  data.  At  high  values  of  anharmonicity  (~9),  the  linear  model  was  able  to  match 

the  universal  cold-curve  (UCC)  fit  of  Rose  et  al  [10]  with  a  value  of  ^proaching  unity. 
The  power  law  does  equally  well  here,  since  for  an  exponent  approaching  unity,  the  linear  and 
power  forms  are  identical.  However,  for  low  anharmonicities  (~3),  no  single  value  of  F^tir'  could 

produce  a  match  to  the  UCC  model  [15,16],  but  instead  required  a  slope  that  was  initially  around 
0,8,  reducing  to  0.55  at  flie  larger  expansions.  Figure  3  compares  the  lattice  potential  energy 
predicted  by  the  UCC  model  of  Rose  to  that  using  the  linear-'n/'  approximation.  In  this  graph,  the 
abscissa,  a,  represents  the  nondimensional,  relative  lattice  spacing  employed  by  Rose  et  al.  [10]. 

By  contrast,  the  power  law  \|;(V)  is  able  to  match  the  UCC  model  closely  with  a  single 
value  of  exponent  x  for  a  given  anharmonicity.  The  value  of  the  exponent  x,  which  will  provide 
a  match  to  the  UCC  model,  is  given  approximately  by 


10 
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Figure  3.  Lattice  binding  energy  as  a  function  of  nondimensionalizcd  lattice  spacing  for  a 
material  Dark  line  is  UCC  model  of  Rose  et  al  [10]  Dashed  lines  are  for  model  of 

Segletes  [15,16],  with  linear  ii;  assumption,  for  two  particular  values  of  and  an 

anharmonicity  of  3. 
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(8) 


where  ti  is  the  value  of  anharmonicity,  as  defined  by  Rose  et  al.  [10].  The  magnitude  of 
anharmonicity  provides  a  metric  on  the  deviation  of  an  actual  crystal’s  behavior  from  the 
harmonic  idealization,  in  which  the  binding  force  of  the  atomic  lattice  would  be  proportional  to 
the  atom’s  displacement  from  the  rest  position.  It  is  later  shown  that  Segletes’  model  [15,16] 
approximates  the  harmonic  solid  described  by  Dugdale  and  MacDonald  [5]  when  the  value  of 
exponent  x  is  selected  as  2/3,  in  accordance  with  the  lower  limit  of  eqn.  (8),  as  Ti-^O.  Figure  4 
compares  the  lattice  binding  energy  predicted  by  the  UCC  model  to  that  of  the  current  model, 
when  a  power-law  \|r(V)  is  employed,  the  exponent  being  selected  in  accordance  with  eqn.  (8). 
Though  not  labeled,  the  dark  curve  represents  the  UCC  model,  while  the  other  five,  essentially 
overlapping,  curves  are  for  the  current  model  with  anharmonicities  of  2,  4,  6,  8,  and  10. 

4.  GRUNEISEN  MODELS 

Certain  mathematical  consequences  of  a  power-law  'Vl/'(V),  eqn.  (7),  which  are  subsequently 
used,  are  easily  derived  and  are  simply  stated  as  follows: 


(roV)v.  =  ^  ’ 

(9) 

< 

II 

• 

(10) 

y  -  n-(}\r/V)  =  (x-n)-(y\r/V)  ;  and 

(11) 

\|/‘[\|/'  -  n-Cllf/Vy]  =  (x  -  n)  (x  -  l)-(\|/'/Vy  . 

(12) 

We  now  intend  to  show  the  primary  result  of  this  paper — ^namely,  that  three  historical 
Griineisen  formulations  (those  of  Slater  [4],  Dugdale  and  MacDonald  [5],  and  the  free-volume 
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Figure  4.  Lattice  binding  energy  as  a  function  of  nondimensionalized  lattice  spacing  for  a 
material.  Dark  line  is  UCC  model  of  Rose  et  al  [10].  Dashed  lines  (essentially  coincident)  are 
for  model  ofSegletes  [15,16],  with  power-law  ii/CV)  assumption,  for  anharmonicity  values  of  2, 
4,  6,  8,  and  10. 
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theory  [6])  may  be  recovered  through  manipulation  of  the  current  model,  if  the  material  is 
idealized  and  higher-order  terms  are  neglected.  Start  with  eqn.  (4)  and  multiply  by  a  power  of 
specific  volume,  V",  where  n  is,  at  this  point,  still  an  arbitrary  constant: 


(0/©of  ln(0/0o)  . 


(13) 


The  first  and  second  derivatives  of  this  equation  may  be  obtained  in  a  tedious  but  straightforward 
manner  as: 


,  E.K^V" 


{[^:+\|/-n(\ir/y)]  ln(0/0o)  +  l}  (14) 


and 


„  E.K^V" 
(P,V'‘)"=  ' 


(0/0o)^  ([A'+2\|/-n(\|f/V)]([is:+\|;'-n(Y/y)]ln(0/0o)  +  l}  + 

■  .  (15) 

[ST + V  -  nUyflV)  ]  -  -  «{Y/V)'  ]  ta(8/0(, )) 


Take  eqn.  (15)  multiplied  by  (-\|r)  and  divide  by  eqn.  (14)  to  obtain: 


^Pc'^")"  ,,  J  n  ,  {[^+V-w(Y/V)f +\ir[V'-n(\|r/V)']}ln(0/0o) 

-vf _ f _  =  2K  +  3vr  -  2n(vr/V)  - - —  .  (16) 


-\|f  _  =  2K  +  3vr'  -  2n(\|/'/V)  - 

(pV”)' 


1  +  [^+\/-n(\|;/V)]ln(0/0o) 


Eqns.  (10^12)  may  be  substituted  into  this  expression  to  eliminate  ti/,  the  left  hand  Y  naay  be 
broken  out  as  V/F,  and  the  equation  may  be  algebraically  manipulated  to  give 
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3x-2n  _  V  (pV’^f  ^  ^  ^ 

'  2  (p,v”y  ■ 

^  ^  _  1  {[^ + (-^  -  »)(V/V)  f  +  ix-  n)(x  -  l)(v/V)^}  ln(0/0(^ 
2  1  +  [i^  +  (A:-n)(\|f/V)]ln(0/0o) 


The  large  term  in  parentheses  at  the  end  of  this  equation  takes  on  the  value  K  at  ambient 
conditions  and  diminishes  with  compression  as  the  logarithm  terms  exert  their  influence.  The 
material  parameter  K,  for  materials  studied  by  Segletes  [15,16],  has  generally  been  close  to  the 
idealized  value  of  unity  [whereupon  (Co/ro)^=Efc].  For  the  sake  of  discussion,  therefore,  consider 
the  large  term  in  parentheses  as  unity,  as  an  approximation  to  the  situation  near  ambient  volume 
[ln(0/0o)«O]  for  an  idealized  material  (iir«l).  Eqn.  (17)  becomes 


3x  -  2n  _  V  (PcV")" 
2 


In  this  equation,  x,  which  denotes  the  power  of  \|r(V),  is  a  material  parameter.  If  we  select  the 
value  of  parameter  n,  which  has  been  arbitrary  to  this  point,  to  be  related  to  x  in  the  following 
way: 


n  ^  3x- 413  ,  (19) 

then  X  may  be  eliminated  from  eqn.  (18)  to  yield 


4-3n  _V  (Pvr 
6  ■  2  (p^V^y 


(20) 


Eqn.  (20)  is  precisely  the  generalized  result  presented  by  Vaschenko  and  Zubarev  [6]  for  the 
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Griineisen  function.  In  this  equation,  if  n=0  {i.e.,  when  x==4/9),  Slater’s  [4]  result  is  obtained. 
If  n=2/3  (i.€.,  when  jc-2/3),  then  the  result  of  Dugdale  and  MacDonald  [5]  is  recovered.  Finally, 
for  the  case  where  n=4/3  when  ;c=8/9),  the  free-volume  theory  [6]  expression  follows.  What 
this  result  tells  us  is  that  any  one  of  the  three  historical  Griineisen  models  may  be  nominally 
applicable,  depending  on  the  properties  of  the  material  of  interest.  In  the  case  of  those  metals 

already  studied  by  Segletes  [15,16],  values  for  ,  which  correspond  to  exponent  x  at  ambient 

conditions  [via  eqn.  (9)],  fell  as  low  as  0.43  («4/9)  for  molybdenum  and  as  high  as  0.88  («8/9) 
for  rubidium,  with  the  vast  majority  of  materials  in  the  0.75-0.85  range. 

5.  NONIDEAL  EFFECTS 

Let  us  now  consider  the  practical  influence  of  the  large  term  in  parentheses  in  eqn.  (17), 
idealized  as  unity  in  the  previous  section.  The  leading  K  term  in  the  braces  affects  the  ambient 
value  of  the  Griineisen  parameter,  while  the  remaining  term  involving  logarithms  governs  the  rate 
at  which  the  Gruneisen  function  changes  with  volume.  If  Segletes’  equation  of  state  is  accepted, 
then  the  fact  that  the  large  term  in  parentheses  diminishes  with  volume  indicates  that  all  of  the 
historical  Gruneisen  formulations  (Slater,  Dugdale-MacDonald,  and  ffee-volume  theory)  will  tend 
to  overestimate  the  rate  of  change  of  the  Gruneisen  function  with  compression  when  using  a 
cold-compression  curve  as  the  baseline. 

Whereas  the  historical  Gruneisen  models  are  based  on  various  theories  of  lattice  vibration, 
the  model  of  Segletes  [15,16],  in  contrast,  is  shown  to  fit  actual  cold  and  Hugoniot  data  well  to 
several  megabars  of  pressure.  Thus,  if  one  accepts  the  validity  of  the  Gruneisen  assumption 
[r=/(V)],  then,  over  the  compression  range  of  actual  thermodynamic  data  to  which  Segletes’ 
model  fits  well,  one  must  accept  as  more  accurate  the  mildly  varying  T  of  that  model  to  the  more 
widely  varying  estimates  produced  by  the  theories  of  Slater,  Dugdale-MacDonald  and/or  the  ffee- 
volume  theory. 

To  see  how  the  models  compare.  Figures  5-7  depict  the  Gruneisen  functional  behavior 
for  a  hypothetical  material  with  parameter  ^=1  and  ro=2.  These  three  figures  depict  four  curves 
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Figure  5.  Griineisen  function  versus  relative  volume  for  hypothetical  material  with  =2  and 
K=l,  using  power-law  with  x=4/9.  Predictions  using  models  of  Slater  [4],  Dugdale  and 
MacDonald  [5],  and  the  free-volume  theory  [6]  are  shown  for  comparison.  Note  that  ambient 
value  associated  with  power-law  matches  that  of  Slater. 


V/Vo 

Figure  6.  Griineisen  function  versus  relative  volume  for  hypothetical  material  with  F^  =2  and 
K=l,  using  power-law  ^^(V),  with  x=2l3.  Predictions  using  models  of  Slater  [4],  Dugdale  and 
MacDonald  [5],  and  the  free-volume  theory  [6]  are  shown  for  comparison.  Note  that  ambient 
value  associated  with  power-law  •\^(V)  matches  that  of  Dugdale  and  MacDonald. 
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Fi^re  7.  Gruneisen  function  versus  relative  volume  for  hypothetical  material  with  Tg  =2  and 
K=l,  using  power-law  iirfyj,  with  x=8/9.  Predictions  using  models  of  Slater  [4],  Dugdale  and 
MacDonald  [5],  and  the  free-volume  theory  [6]  are  shown  for  comparison.  Note  that  ambient 
value  associated  with  power-law  '>^(V)  matches  that  of  the  free-volume  theory. 
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each.  In  each  graph,  there  is  a  Griineisen  curve  corresponding  to  the  power-law  \|/(V),  with 
exponent  x  equal  to  4/9,  2/3,  and  8/9  respectively.  The  resulting  cold  curve  from  the  current 
model  is  used  to  show  what  the  models  of  Slater,  Dugdale-MacDonald,  and  the  free-volume 
theory  would  predict  for  the  behavior  of  the  Griineisen  function,  depicted  in  each  of  the  figures 
by  curves  labeled  S,  D-M,  and  FVT,  respectively.  Per  eqn.  (17),  with  K  equal  to  unity,  the 
ambient  values  for  the  current  thory  match  the  value  of  Slater  when  x=4/9,  that  of  Dugdale  and 
MacDonald  when  x=2l3,  and  that  of  the  free-volume  theory  when  x=S/9. 

To  see  a  comparison  of  the  other  models  to  that  of  Segletes,  for  an  actual  material,  we 
consider  the  case  of  stainless  steel,  because  of  the  availability  of  good  cold-curve  data  out  to 
4  Mbar.  Figure  8  shows  the  cold  curve  and  Hugoniot  for  the  steel.  The  data  points  are  the  same 
as  those  originally  analyzed  by  Segletes  [15,16].  However,  unlike  the  earlier  work,  in  which  a 
linear  \|f(V)  function  was  fit  to  the  data  with  a  slope  of  0.76,  the  cold  curve  is  fit  on  this  occasion 
with  a  power  law,  with  an  exponent,  x,  of  0.72.  The  Hugoniot  fit  is  excellent  as  far  as  data  are 
available  (out  to  2  Mbar).  If  the  Grfineisen  assumption  is  accepted,  then  a  power  law  exponent 
of  0.72  must  provide  an  accurate  reflection  of  the  Griineisen  function  down  to  a  relative  volume 
of  0.68  (where  Hugoniot  data  stop),  and  most  likely  provides  an  accurate  reflection  of  the 
function  down  to  a  relative  volume  of  0.55  (where  cold  data  stop).  Figure  9  shows  the  Griineisen 
function  corresponding  to  this  fit  with  the  solid  curve.  Overlaid  on  the  figure  are  the  results  of 
the  models  of  Slater,  Dugdale  and  MacDonald,  and  the  free-volume  theory,  which  are  derived 
from  the  curvature  of  the  stainless  steel  cold-curve  fit.  As  pointed  out  previously,  the  other 
theories  all  show  a  Griineisen  function  that  decreases  more  rapidly  in  compression  than  the 
current  model  and  data  suggest. 

Of  course,  at  very  high  compressions  and  temperatures,  electronic  and  non-Griineisen 
thermal  effects  will  come  into  play,  and  so  even  if  a  simple  power  law  for  \i/'(V)  might  blend  into 
the  TFD  cold  curve,  it  will  likely  fail  to  capture  the  high-temperature  thermal  effects  properly. 
For  example,  the  shock  Hugoniot  of  most  materials  asymptotes  at  a  V/Vq  value  of  around  0.2. 
Such  an  asymptote  corresponds  to  a  relative  compression,  p,  of  4.  Segletes  [9]  proved  that  the 
asymptotic  value  of  the  Gruneisen  parameter  on  the  Hugoniot  cannot  exceed  (and,  in  fact,  must 
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Figure  8.  Pressure-volume  curves  for  stainless  steel.  Cold  data  to  4  Mbar  shown  with  filled 
symbols.  Hugoniot  data  to  2  Mbar  shown  with  open  symbols.  Curves  represent  cold- 
compression  and  shock-Hugoniot  curves  with  T=1.81  and  power-law  exponent  x,  equal  to 
0.72. 


Figure  9.  Griineisen  function  versus  relative  volume  for  stainless  steel  with  Tq—1.81,  using 
power-law  y^(V),  with  x=-0.72.  Predictions  using  models  of  Slater  [4],  Dugdale  and 
MacDonald  [5],  and  the  free-volume  theory  [6]  are  shown  for  comparison. 
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take  on)  a  value  of  2/p^„,  which  in  this  case,  gives  a  value  of  0.5.  If  the  ambient  value  of  the 
Griineisen  parameter  were  a  typical  value  of  ro=2,  then  this  asymptotic  shock  limit  would  possess 
a  \i/'/\|;o  value  of  0.8  at  a  V/Vq  value  of  0.2.  A  quick  look  at  the  power-law  curves  of  Figure  2 
reveals  tiiat  power-law  \|r(V)  functions  with  exponents  between  0.4  and  1.0  will  not  reach  the 
point  (0.2,  0.8).  Thus,  it  would  seem  that  the  only  way  the  Hugoniot  can  get  to  the  point  (0.2, 
0.8)  while  the  cold  curve  asymptotes  to  the  TFD  curves  is  if  \|r  is  temperature  dependent. 

There  is  some  reason  to  suspect  such  temperature-dependent  behavior,  by  looking,  for 
example,  at  the  aluminum  Hugoniot  fit  of  Segletes  [15,16].  When  using  the  known  ambient 
Grflneisen  value  (Segletes  used  2.39)  with  the  current  model,  the  Hugoniot  prediction  was  quite 
poor.  But  by  using  an  ambient  value  of  half  that  amount,  which  doubles  the  value  of  \|f,  the 
aluminum  Hugoniot  and  cold  curve  were  fit  with  precision  to  9  Mbar  (V/Vq  =  0-4).  One  might 
be  led  to  conclude  that  is,  indeed,  temperature  dependent  and  that,  in  the  case  of  aluminum, 
a  dependence  occiurs  at  relatively  low  temperatures.  Such  a  thermal  dependency  would  be  able 
to  explain  how  the  Hugoniot  can  achieve  relatively  large  asymptotic  values  of  nr/tlto,  while,  at  the 
same  time,  the  cold-curve  value  diminishes  below  this  level.  A  significant  thermal  dependency 
of  the  Griineisen  function  for  aluminum  has  been  reported  by  Grod2ka  [1].  In  particular,  he  cites 
experimental  work  with  shocked  aluminum  foam  that  shows  the  ambient-volume  aluminum 
Griineisen  parameter  to  decrease  from  a  cold  value  of  2.5  to  a  value  below  1.4  at  6500  J/g  of 
specific  internal  energy.  Using  a  value  of  3R/W  to  approximate  the  specific  heat  over  this  range 
(where  R  is  the  universal  constant,  and  W  is  the  atomic  weight  of  aluminum),  6500  J/g 
corresponds  to  a  temperature  of  approximately  7000  K,  and  would  result  from  a  moderate- 
strength  shock  in  the  1.5  Mbar  range. 

As  a  final  note,  another  way  in  which  Segletes’  model  and  those  of  Slater,  Dugdale  and 
MacDonald,  and  the  free-volume  theory  may  be  contrasted  is  accomplished  by  rearranging 
eqn.  (20)  as 
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+  =  0  . 
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As  before,  the  value  of  n  determines  which  of  the  three  older  models  is  being  discussed.  This 
equation  is  first-order  in  By  comparison,  the  current  model,  with  a  power-law  tirCV) 

function,  satisfies  the  following  identity: 

K^(p,V^)  +  [2K  +  x(y\r/V)]'^ipy'‘y  +  =  0  ,  (22) 

which  is  second-order  in  p^^'  These  two  equations  are  obviously  different,  in  that  eqn.  (21) 
contains  only  derivatives  of  while  eqn,  (22)  involves  both  p^^  and  its  derivatives. 

However,  there  is  also  a  great  deal  of  similarity  in  that  not  only  are  both  equations  hyperbolic, 
but  for  the  idealized,  harmonic  case  of  J^T  =  1  and  x  =  n  =  2/3,  eqns.  (21)  and  (22)  differ  only  by 
the  first  term  of  eqn.  (22). 

If  eqn,  (19)  is  used  to  convert  eqn.  (21)  from  n-based  to  x-based  (to  consider  cases  when 
n  ^  x),  it  can  be  shown  (for  K=l)  that  the  transformed  eqn.  (21)  also  differs  from  eqn.  (22)  only 
in  the  first  term  of  eqn.  (22),  regardless  of  the  value  of  x.  It  is  for  for  this  reason  that  at  ambient 
conditions,  when  p^  [and  thus  the  first  term  of  eqn.  (22)]  is  zero,  the  value  of  Fq  from  Segletes 
model  matches  the  older  theories  as  a  function  of  x,  respectively.  We  can  thus  think  of  the  older 
theories  as  each  being  applicable  for  a  particular  class  of  material  (defined  by  material  parameter 
x),  but  at  the  same  time  neglecting  the  p^V^  term  necessary  to  capture  the  empirical  data 
correctly. 

6.  CONCLUSIONS 

In  this  report,  the  recently  proposed  equation  of  state  of  Segletes  [15,16]  has  been 
examined  for  compatibility  with  prior  theory.  In  the  original  work,  an  excellent  match  to 
available  data  and  first-principles  computations  was  achieved  with  the  use  of  an  assumed-linear 
expression  for  the  \|f(V)  function,  defined  as  V/T.  In  the  current  work,  the  use  of  a  power  law 
for  ‘Vir(V)  was  examined  in  lieu  of  the  linear  relationship. 
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A  number  of  positive  outcomes  resulted  with  the  use  of  this  power  law.  First,  the  low- 
volume  limit  for  \|/'(V)  approaches  zero,  which  was  shown  to  be  desirable  in  order  to  have  the 
current  theory  blend  smoothly  into  Thomas-Fermi-Dirac  cold-compression  theory.  Secondly,  over 
the  range  of  available  cold-compression  data,  the  linear  and  power-law  fits  can  be  made  virtually 
indistinguishable,  both  providing  excellent  fits  to  the  data.  Finally,  the  binding  energy  function 
using  a  power-law  y^(V)  assumption  was  able  to  reproduce  the  fit  of  Rose  et  al.  [10]  better, 
especially  in  the  area  of  lattice  expansion.  In  this  regard,  an  expression  given  by  eqn.  (8)  shows 
a  simple  relation  between  material  anharmonicity  and  the  fitted  parameter  of  the  current  model, 
in  order  to  match  the  results  of  Rose  et  al.  closely. 

The  primary  result  of  the  paper  has  been  in  showing  how  Segletes’  model  with  a  power- 
law  iirCV),  for  the  case  of  an  idealized  material,  ignoring  higher-order  terms,  is  able  to  reproduce 
three  historical  Griineisen  function  models,  depending  on  the  value  of  the  fitted  material 
parameter  in  Segletes’  model.  In  particular,  when  the  model  parameter,  which  relates  to  the 
decrease  of  the  Griineisen  function  with  compression,  takes  on  a  value  of  4/9,  the  result  of 
Slater  [4]  is  obtained.  When  the  parameter  is  taken  as  2/3,  the  model  of  Dugdale  and 
MacDonald  [5]  is  recovered.  Finally,  if  the  model  parameter  takes  the  value  of  8/9,  the  firee- 
volume  theory  [6]  expression  follows. 

The  paper  then  has  shown  how  Segletes’  model  differs  from  the  historical  Gruneisen 
models,  if  the  higher-order  terms  are  retained.  It  becomes  apparent  that  the  historical  models  of 
the  Griineisen-function  predict  a  more  rapid  decrease  of  the  Gruneisen  function  with  compression 
than  does  the  current  model.  Some  speculations  have  been  put  forward  here  on  the  possibility 
of  a  temperature-dependent  tjf  function,  which  would,  by  definition,  violate  the  Gruneisen 
hypothesis.  Not  only  would  this  possibility  help  to  explain  some  seeming-incompatibilities 
between  the  behaviors  of  the  high-compression  cold  curve  and  the  asymptotic  Hugoniot,  but  a 
temperature-dependent  \|/'  could  also  explain  the  otherwise  anomalous  behavior  for  aluminum 
cited  by  Segletes  [15,16]  in  the  relative  volume  range  of  0.85  to  unity.  Finally,  the  older 
Gruneisen  models  and  Segletes’  model  have  been  mathematically  expressed  in  a  different  form, 
which  clearly  shows  the  latter’s  distinction  from  the  others.  In  the  case  of  the  comparison 
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between  an  idealized,  harmonic  material  in  the  current  model  to  the  model  of  Dugdale  and 
MacDonald,  the  distinction  originates  from  a  single  pV^'^  term,  which  represents  a  higher-order 
departure  from  the  earlier  theories. 

There  are  several  future  efforts  in  this  area  that  would  be  of  utility.  From  a  position  of 
theoretical  interest,  it  would  be  useful  to  reconcile  the  current  theory  with  the  very-high-pressure, 
equation-of-state  theories,  in  something  more  than  an  empirical  blend.  Barring  this  reconciliation, 
however,  the  current  theory  has  still  shown  its  applicability  into  the  megabar  regime  with 
excellent  results,  and  should  be  of  utility  to  those  modeling  material  deformation  resulting  from 
impact  and  conventional  explosive  events  in  the  ordnance  regime.  Also,  the  area  of  non- 
Grflneisen  (and/or  phase  change)  effects  seems  an  interesting  and  important  area  for  further  work. 
For  many  materials,  the  Griineisen  approach  seems  perfectly  adequate  over  a  large  range  of 
compressions  and  pressures.  However,  there  do  seem  to  be  materials  (aluminum  has  been 
mentioned  in  this  paper)  for  which  non-GrOneisen  effects  at  moderate  pressures  might  provide 
an  adequate  explanation. 
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IOWA  STATE  UNIVERSITY 
ATTN  J  H  ROSE 
AMES  lA  50011 

1  CALIFORNIA  ST  UNIV 
NORTHRIDGE 

ATTN  JON  H  SHIVELY 
10343  4TH  ST  NW 
ALBUQUERQUE  NM  87114 

2  CALTECH 

ATTN  ANDREW  P  INGERSOLL  MS  170  25 
THOMAS  J  AHRENS  MS  252  21 
1201  E  CALIFORNIA  BLVD 
PASADENA  CA  91125 

1  CALTECH 

ATTN  GLENN  ORTON  MS  169  237 
4800  OAK  GROVE  DR 
PASADENA  CA  91007 

1  DREXEL  UNIVERSITY 
ATTN  PHYSICS  DEPT 
32ND  &  CHESTNUT  ST 
PHILADELPHIA  PA  19104 

1  GEORGIA  ESrSHTUTE  OF  TECHNOLOGY 

COMPUTATIONAL  MODELING  CENTER 
ATTN  S  ATLURI 
ATLANTA  GA  30332-0356 


5  JOHNS  HOPKINS  UNIV 
APPLIED  PHYSICS  LAB 
ATTN  TERRY  R  BETZER 
ALVIN  R  EATON 
RICHARD  H  KEITH 
DALE  K  PACE 
ROGER  L  WEST 
JOHNS  HOPKINS  ROAD 
LAUREL  MD  20723 

1  LOUISIANA  STATE  UNIVERSITY 
ATTN  ROBERT  W  COURTER 
948  WYLIE  DR 
BATON  ROUGE  LA  70808 

1  MIT  DEPT  OF  EARTH  ATMOS  AND 

PLANETARY  SCIENCES 
ATTN  HEIDI  B  HAMMELL  54  316 
CAMBRIDGE  MA  02139 

1  NC  STATE  UNIVERSITY 
ATTN  YASUYUKI  HORIE 
RALEIGH  NC  27695-7908 

1  PENNSYLVANIA  STATE  UNIVERSITY 
ATTN  PHYSICS  DEPT 
UNIVERSITY  PARK  PA  16802 

5  SOUTHWEST  RESEARCH  INSTITUTE 
ATTN  C  ANDERSON 
D  LITTLEFIELD 
S  AMULLIN 
B  COUR  PALAIS 
J  WALKER 
PO  DRAWER  28510 
SAN  ANTONIO  TX  78284 

1  STANFORD  UNIVERSITY 
ATTN  BRIAN  PETERSEN 
RAINS  16L 
704  CAMPUS  DR 
STANFORD  CA  94305 

1  TEXAS  A&M  UNIVERSITY 
PHYSICS  DEPARTMENT 
ATTN  DAN  BRUTON 
COLLEGE  STATION  TX  77843-4242 
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1  UC  DAVIS 

INST  OF  THEORETICAL  DYNAMICS 
ATTN  E  G  PUCKETT 
DAVIS  CA  95616 

1  UC  SANTA  BARBARA 

INSTITUTE  FOR  THEORETICAL  PHYSICS 

ATTN  F  GUINEA 

SANTA  BARBARA  CA  93106 

2  UNTV  OF  ALA  HUNTSVILLE 
AEROPHYSICS  RSCH  CTR 
ATTN  GARY  HOUGH 
DAVID  J  LIQUORNIK 

PO  BOX  999 
HUNTSVILLE  AL  35899 

1  UNIV  OF  ALA  HUNTSVILLE 
CIVIL  ENGRNG  DEPT 
ATTN  WILLIAM  P  SCHONBERG 
HUNTSVILLE  AL  35899 

1  UNTV  OF  ALABAMA 

ATTN  A  EUGENE  CARDEN 
86  EASTERN  HILLS 
COTTONDALE  AL  35453 

1  UNIVERSITY  OF  ARIZONA 
ATTN  DAVID  H  LEVY 
TUCSON  AZ  85721 

1  UNIVERSITY  OF  CA  BERKLEY 
MECHNL  ENGRNING  DEPT 
GRADUATE  OFFICE 
ATTNKEZHUNLI 
BERKELEY  CA  94720 

1  UNIVERSITY  OF  CHICAGO 

DEPT  OF  THE  GEOPHYSICAL  SCIENCES 
ATTN  G  H  MILLER 
5734  S  ELLIS  AVE 
CHICAGO  BL  60637 

1  UNIVERSITY  OF  COLORADO 
ATTN  T  MACLAY 
CAMPUS  BOX  431  NNT  341 
BOULDER  CO  80309 


4  UNIVERSITY  OF  DELAWARE 

DEPT  OF  MECHANICAL  ENGINEERING 

ATTN  PROF  J  GILLESPIE 

DEAN  R  B  PIPES 

PROF  J  VINSON 

PROF  D  WILKINS 

NEWARK  DE  19716 

1  UNIVERSITY  OF  MARYLAND 

ATTN  PHYSICS  DEPT  (BLDG  082) 
COLLEGE  PARK  MD  20742 

1  UNIVERSITY  OF  PUERTO  RICO 
DEPT  CHEMICAL  ENGINEERING 
ATTN  L  A  ESTEVEZ 
MAYAGUEZ  PR  00681-5000 

1  UMVERSnY  OF  TEXAS 

DEPT  OF  MECHANICAL  ENGINEERING 
ATTN  ERIC  P  FAHRENTHOLD 
AUSTIN  TX  78712 

1  VIRGINIA  POLYTECHNIC  INSTTIUTE 
COLLEGE  OF  ENGINEERING 
ATTN  R  BATRA 
BLACKSBURG  VA  24061-0219 

1  ADVANCED  TECHNOLOGY  INC 
ATTN  J  ADAMS 

PO  BOX  125 

DAHLGREN  VA  22448-0125 

2  AEROJET  ELECTRO  SYS  CO 
ATTN  WARHEAD  SYSTEMS 
J  CARLEONE 

SKEY 
PO  BOX  296 
AZUSA  CA  91702 

2  AEROJET  ORDNANCE 
ATTN  P  WOLF 
G  PADGETT 
1100  BULLOCH  BLVD 
SOCORRO  NM  87801 
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3  ALLIANT  TECHSYSTEMS  INC 
ATTN  t  HOLMQUIST  MNll  2720 
RSTRYK 

G  R  JOHNSON  MNll  2925 
600  SECOND  ST  NE 
HOPKINS  MN  55343 

1  ALME  AND  ASSOCIATES 
ATTN  MARVIN  L  ALME 
6219  BRIGHT  PLUME 
COLUMBIA  MD  21044-3790 

1  APPLIED  RESEARCH  ASSOC  INC 
ATTN  JEROME  D  YATTEAU 
5941  S  MIDDLEFIELD  RD  SUITE  100 
LITTLETON  CO  80123 

1  APPLIED  RESEARCH  ASSOC  INC 
ATTN  FRANK  MAESTAS 
4300  SAN  MATEO  BLVD  SE 
ALBUQUERQUE  NM  87110 

1  BATTELLE 

ATTN  ROBER  M  DUGAS 

7501  S  MEMORIAL  PKWY  SUITE  101 

HUNTSVILLE  AL  35802-2258 

3  BATTH-LE 

ATTN  S  GOLASKI 
R  JAMESON 
TECHNICAL  LIBRARY 
505  KING  AVE 
COLUMBUS  OH  43201 

3  BOEING  AEROSPACE  CO 

SHOCK  PHYSICS  &  APPLIED  MATH 
ENGINEERING  TECHNOLOGY 
ATTN  R  HELZER 
T  MURRAY 
J SHRADER 
PO  BOX  3999 
SEATTLE  WA  98124 

1  BOEING  HOUSTON  SPACE  STN 
ATTN  RUSSELL  F  GRAVES 
BOX  58747 
HOUSTON  TX  77258 


1  BRIGS  CO 

ATTN  JOSEPH  E  BACKOFEN 
2668  PETERSBOROUGH  ST 
HERNDON  VA  22071-2443 

3  CALIFORNIA  RSCH  &  TECHNOLOGY 
ATTN  R  BROWN 
RFRANZEN 
DORPHAL 
5117  JOHNSON  DR 
PLEASANTON  CA  94566 

1  CALIFORNIA  RSCH  &  TECHNOLOGY 
ATTN  M  MAJERUS 
PO  BOX  2229 
PRINCETON  NJ  08543 

1  COMPUTATIONAL  MECHANICS 
CONSULTANTS 
ATTN  J  A  ZUKAS 
PO  BOX  11314 
BALTIMORE  MD  21239-0314 

1  DEFENSE  TECHNOLOGY  INTERNATIONAL  INC 

ATTN  D  E  AYER 
THE  STARK  HOUSE 
22  CONCORD  ST 
NASHUA  NH  03060 

1  DESKIN  RESEARCH  GROUP  INC 
ATTN  EDWARD  COLLINS 

2270  AGNEW  RD 
SANTA  CLARA  CA  95054 

2  El  DUPONT  DE  NEMOURS  &  CO 
ATTN  L  MINOR 

B  SCOTT 

SECURITY  DIR  LEGAL  DEPARTMENT 
PO  BOX  1635 
WILMINGTON  DE  19899 

1  G  E  DUVALL 

5814  NE  82ND  COURT 
VANCOUVER  WA  98662-5944 
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DYNA  EAST  CORP 
ATTN  P  C  CHOU 
R  aCCARELLI 
WFLIS 

3620  HORIZON  DRIVE 
KING  OF  PRUSSIA  PA  19406 

DYNASEN 

ATTN  JACQUES  CHAREST 
MICHAEL  CHAREST 
MARTIN  LILLY 
20  ARNOLD  PL 
GOLETA  CA  93117 

ELORET  INSTITUTE 
ATTN  DAVm  W  BOGDANOFF  MS  230  2 
NASA  AMES  RESEARCH  CENTER 
MOFFETT  FIELD  CA  94035 

ENIG  ASSOCIATES  INC 
ATTN  J  ENIG 
D  J  PASHNE 
M  COWPERTHWAITE 
SUITE  500 

11120  NEW  HAMPSHIRE  AVE 
SILVER  SPRING  MD  20904-2633 

EXPLOSIVE  TECHNOLOGY 
ATTN  M  L  KNAEBEL 
PO  BOX  KK 
FAIRFIELD  CA  94533 

GB  TECH  LOCKHEED 
ATTN  JAY  LAUGHMAN 
2200  SPACE  PARK  SUITE  400 
HOUSTON  IX  77258 

GB  TECH  LOCKHEED 
ATTN  LUCE.LE  BORREGO  C23C 
JOE  FALCON  JR  C23C 
2400  NASA  ROAD  1 
HOUSTON  TX  77058 

GENERAL  RESEARCH  CORP 
ATTN  A  CHARTERS 
TMENNA 
PO  BOX  6770 

SANTA  BARBARA  CA  93160-6770 


NO.  OF 

COPIES  ORGANIZATION 

2  GRC  INTERNATIONAL 

ATTN  TIMOTHY  M  CUNNINGHAM 
WILLIAM  M  ISBELL 
5383  HOLLISTER  AVE 
SANTA  BARBARA  CA  93111 

5  INST  OF  ADVANCED  TECHNOLOGY 
UNIVERSITY  OF  TX  AUSTIN 
ATTN  STEPHEN  J  BLESS 
JAMES  CAZAMIAS 
HARRY  D  FAIR 
THOMAS  M  KIEHNE 
NICK  LYNCH 
4030-2  W  BRAKER  LN 
AUSTIN  TX  78759 

1  KAMAN  SCIENCES  CORP 
ATTN  DENNIS  L  JONES 
2560  HUNTINGTON  AVE  SUITE  200 
ALEXANDRIA  VA  22303 

8  KAMAN  SCIENCES  CORP 
ATTN  J  ELDER 
RICHARD  P  HENDERSON 
DAVID  A  PYLES 
FRANK  R  SAVAGE 
JAMES  A  SUMMERS 
JAMES  S  WILBECK 
TIMOTHY  W  MOORE 
THY  YEM 

600  BLVD  S  SUITE  208 
HUNTSVILLE  AL  35802 

3  KAMAN  SCIENCES  CORP 
ATTN  SHELDON  JONES 
GARY  L  PADEREWSKI 
ROBERT  G  PONZINI 

1500  GRDN  OF  THE  GODS  RD 
COLORADO  SPRINGS  CO  80907 

4  KAMAN  SCIENCES  CORP 
ATTN  NASIT  ARI 
STEVE  R  DIEHL 
WILLIAM  DOANE 
VERNON  M  SMITH 

PO  BOX  7463 

COLORADO  SPRINGS  CO  80933-7463 
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1  D  R  KENNEDY  &  ASSOC  INC 
ATTN  D  KENNEDY 
PO  BOX  4003 

MOUNTAIN  VIEW  CA  94040 

1  KERLEY  PUBLISHING  SERVICES 
ATTN  G  I  KERLEY 

PO  BOX  13835 

ALBUQUERQUE  NM  87192-3835 

2  KTECH  CORPORATION 
ATTN  FRANK  W  DAVIES 
LARRY  M  LEE 

901  PENNSYLVANIA  NE 
ALBUQUERQUE  NM  87110 

1  LIVERMORE  SOFTWARE  TECH  CORP 
ATTN  J  O  HALLQUIST 
2876  WAVERLY  WAY 
LIVERMORE  CA  94550 

1  LOCKHEED  ENGRNG  &  SPC  SCIENCES 
ATTN  E  CYKOWSKI  MS  B22 
2400  NASA  RD  1 
HOUSTON  TX  77058 

1  LOCKHEED  MARTIN  MISSLE  &  SPACE 
ATTN  WILLIAM  R  EBERLE 
PO  BOX  070017 
HUNTSVILLE  AL  35807 

3  LOCKHEED  MARTIN  MISSILE  &  SPACE 
ATTN  M  A  LEVIN  ORG  81  06  BLDG  598 
M  R  MCHENRY 

T  A  NGO  ORG  81  10  BLDG  157 
111  LOCKHEED  WAY 
SUNNYVALE  CA  94088 

4  LOCKHEED  MISSILE  &  SPACE  CO 
ATTN  JOHN  R  ANDERSON 
WILLIAM  C  KNUDSON 

S  KUSUMI  0  81  11  BLDG  157 
J  PHILLIPS  0  54  50 
PO  BOX  3504 
SUNNYVALE  CA  94088 


1  LOCKHEED  MISSILE  &  SPACE  CO 
ATTN  R  HOFFMAN 
SANTA  CRUZ  FACILITY 
EMPIRE  GRADE  RD 
SANTA  CRUZ  CA  95060 

1  LOCKHEED  NASA  JSC 

SPACE  SCIENCE  BRANCH 
ATTN  JAMES  HYDE 
BOX  58561  MC  B22 
HOUSTON  TX  77258 

1  LOCKHEED  MARTIN  AEROSPACE 
ATTN  D  R  BRAGG 
PO  BOX  5837  MP  109 
ORLANDO  FL  32855 

1  MCDONNELL  DOUGLAS 
ASTRONAUTICS  CO 
ATTN  B  L  COOPER 
5301  BOLSA  AVE 
HUNTINGTON  BEACH  CA  92647 

1  ORLANDO  TECHNOLOGY  INC 
ATTN  DANIEL  A  MATUSKA 
PO  BOX  855 
SHALIMAR  FL  32579 

1  PHYSICAL  SCIENCES  INC 
ATTN  PETER  NEBOLSINE 
20  NEW  ENGLAND  BUS  CTR 
ANDOWER  MA  01810 

3  PHYSICS  INTERNATIONAL 
ATTN  R  FUNSTON 
G  FRAZIER 
L GARNETT 
PO  BOX  5010 
SAN  LEANDRO  CA  94577 

1  RAYTHEON  ELECTRONIC  SYSTEMS 
ATTNRKARPP 
50  APPLE  HILL  DRIVE 
TEWKSBURY  MA  01876 

1  ROCKWELL  INTERNATIONAL 
ROCKETDYNE  DIVISION 
ATTN  J  MOLDENHAUER 
6633  CANOGA  AVE  HB  23 
CANOGA  PARK  CA  91303 
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1  ROCKWELL  INTERNATIONAL 
ROCKETDYNE  DIVISION 
ATTN  H  LEEER 
16557  PARK  LN  CIRCLE 
LOS  ANGELES  CA  90049 

1  ROCKWELL  MISSEE  SYS  DIV 
ATTN  T  NEUHART 
1800  SATELLITE  BLVD 
DULUTH  GA  30136 

1  S  CUBED 

ATTN  R  SEDGWICK 
PO  BOX  1620 
LA  JOLLA  CA  92038-1620 

2  SAIC 

ATTN  JAMES  FURLONG 
GREGORY  J  STRAUCH 
1710  GOODRIDGE  DR 
MCLEAN  VA  22102 

1  SAIC 

ATTN  MICHAEL  W  MCKAY 
10260  CAMPUS  POINT  DR 
SAN  DEGO  CA  92121 

1  SHOCK  TRANSENTS  INC 
ATTN  DAVID  DAVISON 
BOX  5357 

HOPKINS  MN  55343 

2  SIMULATION  &  ENG  CO  INC 
ATTN  ELSA  I  MULLINS 
STEVEN  E  MULLINS 

8840  HWY  20  SUITE  200  N 
MADISON  AL  35758 

2  SOUTHERN  RESEARCH  INSTirUTE 
ATTN  LINDSEY  A  DECKARD 
DONALD  P  SEGERS 

PO  BOX  55305 

BIRMINGHAM  AL  35255-5305 

3  SRI  INTERNATIONAL 
ATTN  JAMES  D  COLTON 
DR  L  SEAMAN 
DONALD  A  SHOCKLEY 
333  RAVENSWOOD  AVE 
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MENLO  PARK  CA  94025 

2  TELEDYNE  BROWN  ENGR 
ATTN  JIM  W  BOOTH 
MARTIN  B  RICHARDSON 
PO  BOX  070007  MS  50 
HUNTSVELE  AL  35807-7007 

1  ZERNOW  TECHNICAL  SVCS  INC 
ATTN  LOUIS  ZERNOW 
425  W  BONITA  AVE  SUITE  208 
SAN  DIMAS  CA  91773 

ABERDEEN  PROVING  GROUND 

26  DE,  USARL 

ATTN:  AMSRL-WM,  I  MAY 

AMSRL-WM-T,  W  F  MORRISON 
AMSRL-WM-TA, 

W  GEUCH 
W  BRUCHEY 
JDEHN 
GFEBEY 
W  A  GOOCH 
HW  MEYER 
E  J  RAPACKI 
AMSRL-WM-TB, 

RFREY 

J  STARKENBERG 
AMSRL-WM-TC, 

W  S  DE  ROSSET 
TWBERKE 
R  COATES 
F GRACE 
KKIMSEY 
M  LAMPSON 
DSCHEFFLER 
S  SCHRAML 
B  SORENSEN 
R  SUMMERS 
W  WALTERS 

AMSRL-WM-PB,  A  ZIELINSKI 
AMSRL-WM-PD,  G  GAZONAS 
AMSRL-WM-WD,  A  PRAKASH 
AMSRL-WM-MF,  A  RAJENDRAN 
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13  DIR,  USARL 

ATTN:  AMSRL-WM-TD, 

A  M  DIETRICH 
KFRANK 
J  HARRISON 
MRAFTENBERG 
G  RANDERS-PEHRSON 
J  SANTIAGO 
M  SCHEIDLER 
S  SCHOENFELD 
S  SEGLETES  (3  CP) 

J  WALTER 
T  WRIGHT 


NO.  OF 

COPffiS  ORGANIZATION 


NO.  OF 

COPIES  ORGANIZATION 


1  ABTEILUNG  FUER  PHYSIKALISCHE 
CHEMIE 

MONTANUNIVERSrrAET 
ATTN  E  KOENIGSBERGER 
A  8700  LEOBEN 
AUSTRIA 

2  AWE 

ATTN  MICHAEL  GERMAN 
WAYNE  HARRISON 
FOULNESS  ESSEX  SS3  9XE 
UNITED  KINGDOM 

1  BATTELLE  INGENIEUTECHNIK  GMBH 
DUESSELDORFFER  STR  9 
ESCHBORN  D  65760 
GERMANY 

1  BHABHA  ATOMIC  RESEARCH  CENTRE 
HIGH  PRESSURE  PHYSICS  DIVISION 
ATTN  N  SURESH 
TROMBAY  BOMBAY  400  085 
INDIA 

1  BULGARIAN  ACADEMY  OF  SCIENCES 
SPACE  RESEARCH  INSTITUTE 
ATTN  VALENTIN  GOSPODINOV 
1000  SOFIA  PO  BOX  799 
BULGARIA 

1  CANADIAN  ARSENALS  LTD 
ATTN  P  PELLETIER 
5  MONTEE  DES  ARSENAUX 
VBLLIE  DE  GRADEUR  PQ  J5Z2 
CANADA 

1  CEA 

ATTN  ROGER  CHERET 
CEDEX  15 

313  33  RUE  DE  LA  FEDERATION 

PARIS  75752 

FRANCE 

1  CEA 

aSI  BRANCH 
ATTN  PATRICK  DAVID 
CENTRE  DE  SACLAY  BP  28 
GIF  SUR  YVETTE  91192 
FRANCE 


1  CEA/CESTA 

ATTN  ALAIN  GEILLE 
BOX2LEBARP  33114 
FRANCE 

3  CENTRE  D’ETUDES  DE  GRAMAT 
ATTN  SOLVE  GERARD 
CHRISTIAN  LOUPIAS 
PASCALE  OUTREBON 
GRAMAT  46500 
FRANCE 

2  CENTRE  D’ETUDES  DE  LIMEIL-VALENTON 
ATTN  CHRISTIAN  AUSSOURD 
JEAN-CLAUDE  BOZIER 

SAINT  GEORGES  CEDEX 
VTLLENEUVE  94195 
FRANCE 

3  CENTRE  D’ETUDES  DE  VAUJOURS 
ATTN  PLOTARD  JEAN-PAUL 
ERIC  BOTTET 

TAT  SIHN  VONG 
BOITE  POSTALE  NO  7 
COUNTRY  77181 
FRANCE 

1  CENTURY  DYNAMICS  LTD 
ATTN  NIGEL  FRANCIS 
12  CITY  BUSINESS  CTR 
BRIGHTON  RD  HORSHAM 
WEST  SUSSEX  RH1354BA 
UNITED  KINGDOM 

1  CRC  RESEARCH  INST  INC 

ATTN  MASAHIDE  KATAYAMA 
STRUCTURAL  ENGINEERING  DEPT 
1  3  316  NAKASE  MIHAMA  KU 
CHIBA  SHI  261  01 
JAPAN 

1  DAT  ETBS  CETAM 

ATTN  CLAUDE  ALTMAYER 
ROUTE  DE  GUERRY  BOURGES 
18015 
FRANCE 
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1  DEFENCE  RESEARCH  AGENCY 
ATTN  W  A  J  CARSON 
CHOBHAM  LANE  CHERTSEY 
SURREY  KT16  OEE 
UNTIED  KINGDOM 

1  DEFENCE  RSCH  ESTAB  SUFFIELD 
ATTN  D  MACKAY 

RALSTON  ALBERTA  TOJ  2NO  RALSTON 
CANADA 

1  DEFENCE  RSCH  ESTAB  SUFFIELD 
ATTN  CHRIS  WEICKERT 
BOX  4000  MEDICINE  HAT 
ALBERTA  TIA  8K6 
CANADA 

1  DEFENCE  RSCH  ESTAB  VALCARTIER 
ATTN  N  GASS 
PO  BOX  8800 

COURCELETTE  PQ  GOA  IRO 
CANADA 

4  DEPARTMENTO  DE  QUIMICA  FISICA 
FACULTAD  DE  CBENCIAS  QUIMICAS 
UNIVERSIDAD  COMPLUTENSE  DE 
MADRID 

ATTN  V  G  BAONZA 

M  TARAVILLO 

MCACERAS 

JNUNEZ 

28040  MADRID 

SPAIN 

1  DYNAMIC  RESEARCH  AB 
ATTN  AKE  PERSSON 
PARADISGRAND  7 
SODERTALJE  S  151  36 
SWEDEN 

4  ERNST  MACH  INSTTrUT 
ATTN  VOLKER  HOHLER 
E  SCHMOLINSKE 
E  SCHNEIDER 
KTHOMA 
ECKERSTRASSE  4 
D-7800  FREIBURG  1  BR  791  4 
GERMANY 


1  ESTEC  CS 

ATTN  DOUGLAS  CASWELL 
BOX  200  NOORDWUK 
2200  AG 
NETHERLANDS 

1  EUROPEAN  SPACE  AGENCY 
ATTN  WALTER  FLURY 
ESOC5 

ROBT  BOSCHT  STRASSE 
DARMSTADT  6100 
GERMANY 

2  EUROPEAN  SPACE  AGENCY  ESTEC 
ATTN  LUCY  BERTHOUD 
MICHEL  LAMBERT 

POSTBUS  BOX  299  NOORDWUK 

NL2200  AG 

NETHERLANDS 

3  FRENCH  GERMAN  RESEARCH  INST 
ATTN  HARTMUTH  F  LEHR 

ROLF  HUNKLER 
ERICH  WOLLMANN 
POSTFACH  1260 
WEIL  AM  RHEIN  D-79574 
GERMANY 

3  FRENCH  GERMAN  RESEARCH  INST 
ATTN  HANS-JURGEN  ERNST 
FRANCIS  JAMET 

PASCALE  LEHMANN 

CEDEX  5  5  RUE  DU  GENERAL  CASSAGNOU 

SAINT  LOUIS  68301 

FRANCE 

4  HIGH  ENERGY  DENSITY  RESEARCH  CTR 
ATTN  VLADIMIR  E  FORTOV 
GENADH  I  KANEL 

V  A  SKVORTSOV 
0  YU  VOJOBIEV 
IZHORSKAJA  STR  13/19 
MOSCOW  127412 
RUSSIAN  REPUBLIC 
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5  lABG 

ATTN  H  J  RAATSCHEN 
WSCHITTKE 
F  F  SCHARPPF 
MBORRMANN 
H  G  DORSCH 
EINSTEINSTRASSE  20 
D  8012  OTTOBRUN  B  MUENCHEN 
GERMANY 

1  INSTITUTE  FOR  PROBLEMS  OF  STRENGTH 

ATTN  G  STEPANOV 
UMIRYAZEVSKAYU  STR  2 
252014  KIEV 
UKRAINE 

1  INSTITUTE  OF  CHEMICAL  PHYSICS 
ATTN  A  YU  DOLGOBORODOV 
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